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Susceptibility tensors in the frequency domain 

The susceptibility tensors in the frequency domain arise when the electric field E a (t) of the light 
is expressed in terms of its Fourier transform E a (u>), by means of the Fourier integral identity 

/oo 
E a (w) exp(-twt) oluj = g-^E^it), (1') 
-oo 

with inverse relation 

E a {uj) = ^J E a (r) exp(^r) dr = $[E a ](w). (1") 

This convention of inclusion of the factor of 2tt, as well as the sign convention, is commonly used 
in quantum mechanics; however, it should be emphasized that this convention is not a commonly 
adopted standard in optics, neither in linear nor in nonlinear optical regimes. 

The sign convention here used leads to wave solutions of the form f(kz-Lut) for monochromatic 
waves propagating in the positive z-direction, which might be somewhat more intuitive than the 
alternative form f(cot — kz), which is obtained if one instead apply the alternative sign convention. 

The convention for the inclusion of 2ir in the Fourier transform in Eq. (1") is here convenient 
for description of electromagnetic wave propagation in the frequency domain (going from the time 
domain description, in terms of the polarization response functions, to the frequency domain, in 
terms of the linear and nonlinear susceptibilities) , since it enables us to omit any multiple of 2tt of 
the Fourier transformed fields. 

First order susceptibility tensor 

By inserting Eq. (1') is inserted into the previously obtained 1 relation for the first order, linear 
polarization density, one obtains 

/oo 
R$,(T)E a (T,t-T)dT, 
-oo 

= {express i£ Q (r, t — r) in frequency domain} 

/OO /'OC 
4aW / E a (r,u))exp[-iu)(t-T)]duidT, 
- oo J — oo 

= {change order of integration} 

/OO /'OO 
/ ^a( T )-^a( r ) a; ) ex P(^ T ) dr exp(— iut) du, 
-oo J — oo 

/oc- 
X^2(-^; w)E a (r, u) exp(-iwi) duj, 
-OO 

where the linear electric dipolar susceptibility, 

p OO 

X$(-w;w) = / R { ^{t) cxp^r) dr - 3ti#)](w), (3) 



(2) 



1 Expressions for the first order, second order, and nth order polarization densities were obtained 
in lecture two. 
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was introduced. In this expression for the susceptibility, ui a — u>, and the reasons for the somewhat 
peculiar notation of arguments of the susceptibility will be explained later on in the context of 
nonlinear susceptibilities. 

Second order susceptibility tensor 

In similar to the linear susceptibility tensor, by inserting Eq. (1') into the previously obtained 
relation for the second order, quadratic polarization density, one obtains 

/OO pOO pOO pOO 
/ / / < 2 ^(n,r 2 )£; Q (r, Wl )^(r, W2 ) 
-oo J — oo J —oo J — OO 

x exp[— i(uji(t — n) + uj 2 (t — r 2 ))] dr\ dr 2 dui duj 2 (4) 

/OO pOO 
/ X^lpi-Ua] wi, Lu 2 )E a {r, wi)£/5(r, lu 2 ) exp[-i {uj x + Lu 2 )t] dui duj 2 
-oo J — oo s v ' 

where the quadratic electric dipolar susceptibility, 



/oo />oo 
/ R ^( T i,T2) exp[i(win + uj 2 t 2 )} dn dT 2l (5) 
-oo J — oo 



was introduced. In this expression for the susceptibility, ui a = u>i + lu 2 , and the reason for the 
notation of arguments should now be somewhat more clear: the first angular frequency argument 
of the susceptibility tensor is simply the sum of all driving angular frequencies of the optical field. 

(2) 

The intrinsic permutation symmetry of i2L^g(ri, t 2 ), the second order polarization response 
function, carries over to the second order susceptibility tensor as well, in the sense that 

i. e. the second order susceptibility is invariant under any of the 2! = 2 pairwisc permutations of 
(a,wi) and (/3,lo 2 ). 

Higher order susceptibility tensors 

In similar to the linear and quadratic susceptibility tensors, by inserting Eq. (1') into the previously 
obtained relation for the nth order polarization density, one obtains 

/oo />oo 
••• / X^ai-aj-w^wi, . . . ,u n )E ai (r,u>i) ■ ■ -E an (r,uj n ) 
-oo J-oc , . 

x exp[— i (uj\ + . . . + u> n ) t] dbJi ■ ■ ■ du> ni 



where the nth order electric dipolar susceptibility, 



(n) 

V 



/OO P oo 

■• / ^ai".a n (Ti } ... ) r n )exp[i(a;iri+. . .-\-uj n r n )] d^ • • • dr ni 
-oo J —oo 

(7) 

was introduced, and where, as previously, 

(jj a = LOi + LU 2 + . . . + LJ n . 

The intrinsic permutation symmetry of i?^a\... Q „ (n, . . . , t„), the nth order polarization re- 
sponse function, also in this general case carries over to the nth order susceptibility tensor as well, 
in the sense that 

( n ) / \ 

Xnaia 2 —a n V ^fi W\ , (J 2 , . . . , U) n ) 

is invariant under any of the n! pairwise permutations of (ai,wi), (a 2 ,u) 2 ), . . ., (a n ,ui n ). 
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Monochromatic fields 

It should at this stage be emphasized that even though the electric field via the Fourier integral 
identity can be seen as a superposition of infinitely many infinitesimally narrow band monochro- 
matic components, the superposition principle of linear optics, which states that the wave equation 
may be independently solved for each frequency component of the light, generally does not hold 
in nonlinear optics. 

For monochromatic light, the electric field can be written as a superposition of a set of distinct 
terms in time domain as 

E(r, t) = Y^ Re [ E ^ exp(-iwfci)], 

k 

with the convention that the involved angular frequencies all are taken as positive, Uk > 0, and 
the electric field in the frequency domain simply becomes a superposition of delta peaks in the 
spectrum, 

1 f°° 

E(r, w) = — / E(r, r) exp(iwr) dr 

27r J-oo 

= {express as monochromatic field} 
1 f°° 

= — ^ / Rc[E Wfc exp(-iw fe T)] exp(io;r) dr 

k "'- oc 
= {by definition} 
1 f°° 

= 4^ J t E "fc ex P(*( w - u k)r) + E* fc exp(i(w + u k )r)] dr 

2 r 1 r°° i r°° 

= o E ^ V~ \ CX P('( W ~ ^fe) T ) dT +E ^ k y~ / exp(i(w + w fe )r) dr 



k _ . ______ . 

=i5(a>— Wfe) =5(_.+w fc ) 

= {definition of the delta function} 

k 

By inserting this form of the electric field (taken in the frequency domain) into the polarization 
density, one obtains the polarization density in the monochromatic form 

P (n) (r,i)= ]T Re[P£>cxpH^t)], 
with complex-valued Cartesian components at angular frequency ui a given as 

( P( Z } )H = 2£ oJ2' • ■£x& ) ia 2 ..-a B (-W (J ;t«;i,W2, . . . , W n )(E Ul ) ai (E^ ) a2 ■ ■ • (-B w „)a„ 

+ X/^ _a 2 ---c_ n ( — W CT ;W2, • • • 7 w n) (^w 2 )ai (^wi )a 2 ' ' ' {^u n )a„ 

+ all other distinguishable terms 

where, as previously, t_v = u>i + _ + • • • + oj n . In the right hand side of Eq. (8), the summation 
is performed over all distinguishblc terms, that is to say over all the possible combinations of 
u>i,u)2, ■ ■ ■ ,co n that give rise to the particular u) a . Within this respect, a certain frequency and 
its negative counterpart are to be considered as distinct frequencies when appearing in the set. 
In general, there are several possible combinations that give rise to a certain uj a ; for example 
(u),co, -co), (lj, — u),ui), and (—u),u,lo) form the set of distinct combinations of optical frequencies 
that give rise to optical Kerr-effect (a field dependent contribution to the polarization density at 

UJ a = LU + LU — LU = Uj). 

A general conclusion of the form of Eq. (8), keeping the intrinsic permutation symmetry in 
mind, is that only one term needs to be written, and the number of times this term appears 
in the expression for the polarization density should consequently be equal to the number of 
distinguishable combinations of u)\, 0J2, ■ ■ ■ , w B . 
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Convention for description of nonlinear optical polarization 

As a "recipe" in theoretical nonlinear optics, Butcher and Cotter provide a very useful convention 
which is well worth to hold on to. For a superposition of monochromatic waves, and by invoking 
the general property of the intrinsic permutation symmetry, the monochromatic form of the nth 
order polarization density can be written as 

{ P ^)n = £0^2 ■ ■■^2^2 K (- U} ^ U U ■ ■ ■,W n )x^a ) 1 -a n (-W ( r;Wl, ■ • ■ , W n ) (E Ul ) ai ■ ■ ■ (E Un ) an . 

(9) 

The first summations in Eq. (9), over ai, . . . ,a n , is simply an explicit way of stating that the 
Einstein convention of summation over repeated indices holds. The summation sign ^2 U , however, 
serves as a reminder that the expression that follows is to be summed over all distinct sets of 
u)\, . . . , u n . Because of the intrinsic permutation symmetry, the frequency arguments appearing in 
Eq. (9) may be written in arbitrary order. 

By "all distinct sets of u)\, . . . , w„", we here mean that the summation is to be performed, as 
for example in the case of optical Kerr-effect, over the single set of nonlinear susceptibilities that 
contribute to a certain angular frequency as (— oj; w, u), — ui) or (—uj;uj,—uj,uj) or (—oj;—oj,uj 7 uj). 
In this example, each of the combinations are considered as distinct, and it is left as an arbitary 
choice which one of these sets that are most convenient to use (this is simply a matter of choosing 
notation, and does not by any means change the description of the interaction) . 

In Eq. (9), the degeneracy factor K is formally described as 

K(-ujy,u Jl ,...,uj n )=2 l+m - n p 

where 

p = the number of distinct permutations of u>i, u>2, ■ ■ ■ , w\, 
n = the order of the nonlinearity, 

m = the number of angular frequencies ojk that are zero, and 

l= (l, if u>„ ± 0, 
[ 0, otherwise. 

In other words, m is the number of DC electric fields present, and I = if the nonlinearity we are 
analyzing gives a static, DC, polarization density, such as in the previously (in the spring model) 
described case of optical rectification in the presence of second harmonic fields (SHG). 

A list of frequently encountered nonlinear phenomena in nonlinear optics, including the degen- 
eracy factors as conforming to the above convention, is given in Butcher and Cotters book, Table 
2.1, on page 26. 

Note on the complex representation of the optical field 

Since the observable electric field of the light, in Butcher and Cotters notation taken as 
E (M) = \ c M-i^kt) + E* fc exp(iwfct)], 

LU k >0 

is a real-valued quantity, it follows that negative frequencies in the complex notation should be 
interpreted as the complex conjugate of the respective field component, or 

E -^ fc = E* fc . 

Example: Optical Kerr-effect 

Assume a monochromatic optical wave (containing forward and/or backward propagating compo- 
nents) polarized in the xy-plane, 

E(z,i) = Re[E w (z) exp(-iw t )] e M 3 , 
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with all spatial variation of the field contained in 

E u (z)=e x EZ(z)+e y Ey(z)&C 3 . 

Optical Kerr-effect is in isotropic media described by the third order susceptibility 

(3) / \ 

X^/3 7 (-w;w,a;,-a;), 

with nonzero components of interest for the zy-polarized beam given in Appendix 3.3 of Butcher 
and Cotters book as 

Y (3) = Y (3) Y (3) = (3) = f intr. perm. symm. 1 (3) (3) (3) (3) 

Axxxx Ayyyyi Axxyy Ayyxx > (a , U>) ^± (0 , U)) J Axyxy Ayxyxi Axyyx Ayxxyi 



with 



Y (3) (3) , (3) , (3) 

Axxxx Axxyy 1 Axyxy 1 Axyyx' 



The degeneracy factor K(—ui;u),u), —w) is calculated as 

K(-u; u, w, -u) = 2 l+m - n p = 2 1+ °- 3 3 = 3/4. 

From this set of nonzero susceptibilities, and using the calculated value of the degeneracy factor in 
the convention of Butcher and Cotter, we hence have the third order electric polarization density 
at oj a — oj given as p(")(r, t) = Rc[p£™' exp(— iu>t)}, with 

= {Using the convention of Butcher and Cotter} 

= 2 ^ £ ° 4 S S 51 x2/3 7 (-W; w > w > -w)(^a,)a(-B«) J 8(^-u;)7 
A* a /3 7 

= {Evaluate the sums over (x,y,z) for field polarized in the plane} 

Q 

= £n -| e [y (3) E X E X E X + y (3) E v E y E x + y (3) E v E x E v + y (3) E x E v E y 1 

4-e [y (3) E v E y E y + y^ 3) E X E X E V + y' 3) E x E v E x + E v E x E x W 

1 J/ VAyyyy u) UJ — 1 Ayxxy ij w -w 1 Ayxyx u uj — lu 1 Ayyxx u u — wJJ 

= {Make use of E_ w = E* and relations = X$L, etc.} 

Q 

— «r n _|e [y (3) /^l/^l 2 + y (3) E y2 E x * + y (3) l/^l 2 /? 21 4- y (3) /^l/^l 2 ! 

x \&xxxx u\ u\ ~ Axyyx-^u ' Ax yxy I ^oj I ~ Axxyy^w \ J -'u I J 

x e [ Y P) K^l 2 + Y (3) E x2 E y * 4- y (3) \E X \ 2 E V 4- y (3) E^li^Pn 

^ ^ylAxxxaj^wl-^wl ~ Axyyx^LO ^u) ' Axyxy \ 1 -'ixj I ~ Axxyy-^LU \ 1 -'ixj \S 

= {Make use of intrinsic permutation symmetry} 

= ^^{^[(xiiiral^J + 2xiiyj/l^l + (Xxxxx ~ ^Xxxyy)^u E u 

e hY (3) l/? y I 2 4- 2y (3) IP^I 2 )^ 4- fy( 3 ) - 2y (3) )E x2 E v * 

^VlXAxxxxl^ujl 1 ^Axxyy I I /-^w 1 VAxxxa; ^Axxyy j-^uj ^lj • 

For the optical field being linearly polarized, say in the x-direction, the expression for the polar- 
ization density is significantly simplified, to yield 

p(f)= £o (3/4)e K xilJ^| 2 ^, 

i. e. taking a form that can be interpreted as an intensity-dependent (~ |^| 2 ) contribution to the 
refractive index (cf. Butcher and Cotter §6.3.1). 
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